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Abstract. In this paper we investigate the properties of nodal structures in 
random wave fields, and in particular we scrutinize their recently proposed 
connection with short-range percolation models. We propose a measure 
which shows the difference between monochromatic random waves, which are 
characterized by long-range correlations, and Gaussian fields with short-range 
correlations, which are naturally assumed to be better modelled by percolation 
theory. We also study the relevance of the quantities which we compute to the 
probability that nodal lines are in the vicinity of a given reference line. 
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1. Introduction 

The nodal domains of a (real) wavefunction are regions of equal sign, and are bounded 
by the nodal lines where the wavefunction vanishes. Even a superficial look at the 
nodal domains of a quantum wavefunction reveals the separable or chaotic nature 
of the quantum system pQ. In separable systems, one observes a grid of intersecting 
nodal lines, and consequently a checkerboard- like nodal domain pattern. In (quantum) 
chaotic systems on the other hand, the nodal domains form a highly disordered 
structure, resembling the geometry found in critical percolation. Blum et al [5] 
argued that also the statistics of the number of nodal domains reflects the fundamental 
difference between separable and chaotic quantum systems. Bogomolny and Schmit 
[3] conjectured that the nodal domain statistics of chaotic wavefunctions in two 
dimensions can be deduced from the theory of critical percolation. They built a 
percolation model for the nodal domains which allowed them to calculate exactly the 
distribution of numbers of domains. Its predictions have been confirmed numerically as 
far as nodal counting and the area distribution of nodal domains are concerned. While 
quantum wavefunctions display long-range correlations, the critical percolation model 
assumes that such correlations can be neglected on distances of the order of a wave 
length. One may thus expect that some nodal properties in real wavefunctions are not 
well described by critical percolation. The main motivation of the present study was to 
investigate the limits of applicability of the short-range percolation model. The object 
we will address is related to the distribution of shapes of nodal lines in the random 
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wave ensemble. To be precise, we will calculate approximately the probability, that a 
nodal line matches a given reference line up to a given precision e. 

We will examine the statistics of nodal lines within the monochromatic random 
wave model, which is a good description of the eigenfunctions of a quantum billiard 
in the semiclassical limit 0]. The monochromatic random wave ensemble consists of 
solutions of the Helmholtz wave equation for a fixed energy E = k 2 

-V 2 $ = /c 2 $. (1) 

Furthermore the random function $ is picked up from a Gaussian distribution, 
i.e. higher order correlations of $ can be expressed through the two-point 
correlation function G\(r) = ($(r)$(r')) by virtue of Wick's theorem. A convenient 
representation of $ is given by the superposition of cylindrical waves with Gaussian 
distributed amplitudes 

$(r, 6) = A m J m (kr) exp{im8) (2) 

m 

where J m {x) are the Bessel functions of the first kind, and r, 9 is the position in polar 
coordinates. The Gaussian random variables obey A* m = (— l) m A m to render $ real, 
and have correlations (A^Am) = S m .m, to, to > 0. Using the addition theorem for the 
Bessel functions, one finds for the two-point correlation function 

Gx(r) = (*(r)*(0> = Mk\r - r>\) ~ ^y -r\-n/A) _ 

Wk\r' — r\ 

It displays in fact long-range correlations, which decay with a power law. In order 
to access the relevance of the long-range correlations, we compare the monochromatic 
random wave ensemble with another Gaussian ensemble of random functions, which, 
however, does not have long-range correlations, and is characterized by the correlation 
function 

G (r) = exp(-fcV/4). (4) 

For the latter ensemble the applicability of the critical (short-range) percolation 
picture is evident jSJ HJ , since the sign of the random function is not significantly 
correlated for distances r ^> Figure ^ shows the spatial correlation functions 

Go, Gi. 

We now briefly introduce the central object of this article. A detailed derivation 
will be given in sectional Consider a smooth, closed reference curve r(s) in the plane, 
which is parametrized by its arclength s. The integral of the square of the amplitude 
of a random function $(r) along this curve 

X = \J d S $(r( S )) 2 (5) 

is itself a random variable. It samples the function not only at a discrete set of points, 
but along a one-dimensional subset of the plane. It should be well suited to detect 
the long-range correlations of the random field $ 

Now assume that $ has a nodal line very close to the given reference line. Then 
X will be small in a sense which will be explained later. Thus, by calculating the 
distribution of X, its cumulants or moments, one obtains the relative importance 
of the given reference line r(s). We will perform these computations for a circular 
reference line both for the random wave ensemble and for the short-range ensemble 
defined above. We will study in particular the scaling properties of the cumulants of X 
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Figure 1. The spatial correlation functions G± (solid), and Go (dotted) as a 
function of kr. 



as functions of the radius (typical size) of the reference curve. We shall show that they 
obey a scaling law which distinguishes clearly between the short-range ensemble and 
the monochromatic random wave ensemble. To understand the significance of these 
cumulants, we shall consider an approximate expression for the probability that a nodal 
line is found inside a strip of width e about the reference line. We shall show that 
this function, when expanded in powers of e _1 generates the X cumulants. Its scaling 
properties with the size parameters, however, are less sensitive to the correlations 
assumed for the underlying random functions model. Strictly speaking, the function 
we compute is a measure of the intensity of fluctuations of the field <& along the line, 
which are certainly small when a nodal line approximates the reference line, but can 
also be small if different nodal lines which just avoid crossing, are within e from the 
reference line. Since near avoided crossings have low probability (see [7j) we believe 
that the function we compute is closely related to the true probability. 

The rest of the paper is organized in the following way. The next section describes 
in detail the new concept which we introduce to the morphological study of nodal 
lines, that is, the density of line shapes. Once this is done, a formal expression for the 
density, expressed in terms of X is provided, and computed explicitly for particular 
shapes — circles (section 01 within a reasonable and calculable approximation. These 
densities are evaluated for random waves, and for the short-range ensemble. 

2. The density of nodal line shapes 

We consider two-dimensional, Gaussian random fields, and a prescribed (closed) 
reference line. We shall propose a proper definition of the density of nodal lines 
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which match the reference line in a random field (or equivalently, the probability 
that a nodal line with a prescribed form shows up in a Guassian random field). 
Compared to problems, where the density of (critical, nodal) points of a Gaussian 
field is calculated [5] ITT] , we enter here a new dimension and consider the 
density of one-dimensional strings instead of zero-dimensional, point-like objects. In 
order to obtain a well-defined and finite theory, we have to regularize the theory by 
dilating the reference curve to a thin tube with constant thickness d and compute the 
probability, that a nodal line is completely inside this tube — see Figure (J2J. Assume 
now, that a function $(r) has a nodal line close to a reference curve r(s), where s 
denotes the arclength. The normal distance rj(s) of the nodal line from the reference 
curve can be obtained via linearization 

$(r + r/n) ps r]d n $(r) + $(r) = (6) 

yielding 

The unit vector n(s) is normal to the curve r(s). d n denotes the corresponding 
normal derivative. The probability, that a nodal line lies in a sharp tube |f?(s)| < d 
is, although well defined, not accessible by analytical means. At this point we must 
resort to further approximations, which will eventually lead us to a tractable model, 
at the cost of losing the rigour of the original object defined above. As a first step we 
replace the box shaped cross section by a smooth Gaussian and consider instead the 
expectation value 

P e = (ew(-±-Jdsr, 2 ^\ (8) 

where J ds is the line integral along the reference line, and e ~ d 3 . However, even the 
computation of this quantity poses unsurmountable difficulties, i] — <p/d n (j) is a ratio 
of two (in general non-independent) Gaussian variables, which is itself non-Gaussian. 
In order to obtain a tractable expression we approximate the integral w = J ds rj 2 by 
a mean-field type expression 

where the latter step requires isotropy of the distribution of the random field <E>. The 
final approximation for the shape probability now reads 

^(^(-^iall^ Kww <10) 



or 



<(V$) 2 )>/2 

where B is an integral operator with (symmetric) kernel 

B(s, s') = (<S>(r(s))<S>(r(s>)) = G(\r(s) - r(s')\) (12) 

and /3 M are the corresponding eigenvalues. The operator B is the correlation function 
of the field <&, restricted to the given curve. It is positive semi-definite and has a finite 
trace J dsB(s, s) — L, thus its eigenvalues /3 M > have an accumulation point at 
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Figure 2. A section of the nodal set of a random wavefunction. One of its 
nodal lines lies within the prescribed thin circular tube, i.e. this configuration 
contributes to the density p. 



zero. The final expression for the logarithm of probability is the starting point of 
our investigation, ft should reflect the relevant features of the inaccessible hard-tube 
probability, and is an interesting object in its own right§. It takes into consideration 
the random field $ along the whole reference curve r(s). We remark here again that 
the final approximation for the probability only tests whether 3> is small along the 
given curve-it is not able to resolve nearly avoided intersections which are placed next 
to the reference curve. 

F(e) is the generating function for the cumulants of the random variable X = 
(1/2) J ds^ 

F( e ) = log(exp(- e ~- 1 X))= -} 7[ (-lY~e- V {X v ) c (13) 

where the expansion parameter is e = e ((V<£ > ) 2 ^ /2. It is also the generating function 
of the traces of powers of the operator B. In fact, expanding F(e) in terms of e, i.e. 
for large e , one finds 

F(e)=-i^log(l + e- 1 /3 m )= £ i-(-l)^-^(/? m )" (14) 

m u— 1,2,3... m 

§ Private communication with J Hannay. 
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Comparing the two expansions, we see that 



v\ 



(*">c = g ( 15 ) 

m 

As mentioned in the introduction our goal is to compare two different Gaussian 
random fields in two dimensions with correlation functions ($(r)$(0)) = G(r) = 
J d 2 p G(p) exp(ip • r), namely 

G 1 (r) = J (kr) 

G (r) = exp(-fcV/4). (16) 

Gi is the correlation function of the monochromatic random wave ensemble with a 
sharply defined energy k 2 . Consequently, it displays long-range correlations. Go is a 
typical short-range ensemble. Note that the G are normalized such that 

G(0) = (<i> 2 ) = J d 2 P G(p) = l 

-V 2 G(0) = <(V$) 2 ) = J d 2 pp 2 G(p) = k 2 . (17) 

This implies an equal nodal line density (| V$|<5(<I>)) for the long, and short-range 
ensemble. 



3. The density of circular nodal lines 

We consider now the approximate probability (|10fl for circles with radius R. The 
kernel of the operator B reads for the monochromatic random wave ensemble with 
correlation function G\ 

B(6-d') = J (2kRsm(^-^-j) (18) 

where 6, 9' are angles describing positions on the circle. Owing to the rotational 
invariance of the problem, the eigenfunctions of B are exp(im0),m = 0,±1,±2, . . .. 
The eigenvalues of the integral operator are therefore 

r-2n 



(3 m =R dd J (2fc J Rsin (9/2)) exp (im6 
Jo 



= 2nR{J m (kR)f (19) 
The eigenvalues for the short-range ensemble read 



f3 m = 27TRexp(~k 2 R 2 /2)I m (k 2 R 2 /2) » ^ exp J ■ (20) 

Figure 13 shows the eigenvalues of B for a circle with radius kR = 100. The spectrum 
for the random waves has strong fluctuations, whereas the spectrum for the short- 
range ensemble is a smooth (almost) Gaussian. It was mentioned before, that the 
trace obeys J2 m @ m ~ ^ f° r both the short- and the long-range ensemble. 
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Figure 3. The spectrum of B as a function of the order m for kR = 100. Shown is 
the random wave case (symbol +), and the short-range case (points are connected 
to a dotted line). 



4. The random wave case 

We calculate now F(e) — logP e and its large-e expansion for large radii R 3> 1/k. In 
the region to < kR for large kR the Bessel functions are well approximated by (see 

EH) 



J m (m sec/3) 



1/2 



cos (mtan/3 — m/3 — 7r/4). 



(21) 



^7r?7itan/3 / 
By setting to sec [3 = kR, we obtain 

J m (kR) « (2/tt) 1 /2 (( kR f _ m 2)-V4 

x cos (fc-R) 2 — "i 2 — marccos(m/(fci?)) — 7r/4^ . (22) 

In the transition region m « fci?, we approximate the Bessel function J m (kR) in terms 
of an Airy function Ai(x) (see \V2\ ) 



JmikR) « ( - 



1/3 



1/3 



(TO - jfci?) 



(23) 



We can combine both asymptotic expansions into a scaling law with a universal scaling 
function f(x) 

\j m m\-{kR)-^f{ m \ k ^f ) (24) 

Figure 01 shows, that the scaling functions f(x) collapse well for three different values 
of kR = 50, 100, 200. Note that for negative arguments the scaling function f(x) 
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Figure 4. The scaling function f(x) for kR = 50 (red +), kR = 100 (green x), 
and for kR = 200 (blue *). 



is strongly fluctuating. In the subsequent applications, f(x) and its powers will be 
integrated over, and for this purpose, f(x) for x < can be considered as a stochastic 
function. f(x) vanishes exponentially for x -*■ +oo, /(0) = (2/9) 1/3 /r(2/3) = 0.44731, 
and f(x) <~ (— a;) -1 / 4 for a; — > — oo. The eigenvalues of the operator scale according to 

^f(->-('(^)) 3 - 

The leading behaviour of the i/th cumulant of X which is proportional to the trace of 
the vth power of B, as a function of the radius R reads 
1u f°° / / , \ \ 2i/ 



(JT") C = $> m r ~ k- v {kR) 1+v ' 3 / di (/ ((t 2 - l)(feii) 2 / 3 )) (26) 



We find to leading order in kR 

( kR v < 2 

(A^) c - jST" x \ kRlog(kR) v = 2 (27) 

[ (fci?)( 1 +")/ 3 i/>2. 

This scaling behaviour is compared with the corresponding quantity for the short- 
range correlations, where 

(X u ) c ~k- u kR (28) 

for all v > 0. Some remarks are in order. The cumulants {X v ) show a typical critical 
behaviour for the random wave case. Below the critical power v* — 2, the large /en- 
sealing does not differ from the short-range case. At the critical power, logarithmic 
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deviations show up, and above v*, the cumulant (X v ) displays an anomalous scaling 
in kR, different from the non-critical, short-range ensemble. Now we return to the 
shape probability 

logP £ « -( kR ) dtlogfl+uj- 1 (kR) 1 / 3 (f ((t 2 - l)(feE) 2 / 3 ))^ (29) 

where u> = ek 3 /(in) is the dimensionlcss width of the tube around the (here circular) 
reference curve. The limit of small lu corresponds to the cumulant (|26[) for v \ as 
far as the scaling behaviour is concerned. Therefore, logP e cannot be considered as a 
'good' quantity to distinguish between the long-range and the short-range cases — for 
both ensembles, log P e ~ kR. There might be anomalous higher order corrections in 
the random wave case, which are not considered here. On the other hand, a large-e 
expansion (which means arbitrarily wide tubes) yields the sequence of cumulants (|26[) 
for integer v which in fact have characteristic scaling properties for v > 2. Figure 
shows a log-log plot of the third cumulant [y = 3) as a function of the radius kR for 
10 < kR < 200. The slope is 1.34063 (standard error = 0.065%), i.e. confirms the 
predicted exponent 4/3. 

5. Conclusion 

In this paper we compared monochromatic Gaussian random waves and a short-range 
ensemble of random fields by investigating the statistics of (1/2) /ds<i> 2 along a 
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given reference curve. The uth order cumulants of this random variable obey non- 
trivial scaling laws with respect to the linear size of the reference curve (here circles 
of radius R) in case of the long-range random waves. The second-order cumulant 
shows logarithmic deviations from the corresponding scaling behaviour of the short- 
range ensemble. The cumulants of order three and higher have non-trivial exponents. 
Namely, these cumulants scale like R^ 1+1/ ^ 3 , whereas the cumulants for the short- 
range ensemble scale ~ R. The probability that a nodal line lies in a circular tube 
of given thickness e, however, turned out to be a less useful candidate to probe the 
long-range properties of the random functions. The logarithm of the shape probability 
for the short- and the long-range case scale in exactly the same manner, which might 
explain the success of the ad hoc model [3] . 
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